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Abstract 

Vector SO (3) gauged 0(4) sigma models on R3 are presented. The 
topological charge supplying the lower bound on the energy and ren- 
dering the soliton stable coincides with the Baryon number of the 
Skyrmion. These solitons have vanishing magnetic monopole flux. To 
exhibit the existence of such solitons, the equations of motion of one of 
these models is integrated numerically. The structure of the conserved 
Baryon current is briefly discussed. 
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1 Introduction 



The problem of gauging a sigma model has been considered in the past, first 
by Fadde'ev|]J and also in the works of Witten0 and Rubakov||] for the 
SkyrmeQ model on R3, which is esentially the 0(4) sigma model. Rubakov[[| 
in particular considers the properties of the soliton in the SU{2) gauged 
Skyrme model. The purpose of the present work is to construct a topologi- 
cally stable soliton in a particular SO (3) gauged Skyrme model. 

The problem of gauging a Sigma model such that the resulting system 
supports topologically stable finite action or energy solutions was recently 
considered in || and |J, the first for the U(l) gauged CP 1 Grassmanian 
model on IR 2 and the second for S0(2n) gauged Grassmanian models on 

The only limitation of these models is that their definitions are restricted 
to even dimensions and the physically important problem of the soliton of 
the Skyrme model[|J, which is the 0(4) Sigma model defined on IR 3 . Indeed 
with this problem in mind, a peculiar gauging of the 0(3) model on IR 2 was 
proposed in |J. Denoting the Sigma model fields a = (0 a , 3 ) with a = 1, 2, 
and <fi a (j) a = 1, the gauging prescription proposed was stated via the minimal 
coupling in terms of the covariant derivative Di<p a = (Di(j) a , Di(j) 3 ) as: 

A0 Q = d t <p a + A i£ a ^, A0 3 = ^0 3 . (1) 

The topological invariants we sought to employ in || were the integrals 
of densities which were total divergences. Since the topological charge den- 
sity of the 0(3) sigma model is not a total divergence but is only locally a 
total divergence 0] , we ignored the 0(3) model gauged with the U(l) field 
according to (|1|) and proceded|| instead to gauge the CP 1 model. Subse- 
quently however the 0(3) model on R 2 , gauged according to ([!]) was shown 
by Schroers[^| to support topologically stable and even self-dual solutions. 

It is our aim here to gauge the 0(4) gauge model @] on IR 3 employing an 
extended version of the minimal coupling prescription given by ([!]), where 
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the U(l) is replaced by the SO (3) gauge group, and to show that the cor- 
responding gauged system supports stable solitons. While we shall restrict 
our considerations here to this 3 dimensional case, it should be noted that all 
our considerations in the present work can be extended to the d dimensional 
case systematically, with the SO(d) gauging of the 0(d+ 1) Skyrme-Sigma 
model||. The starting point in that case would be the generalisation of the 
minimal coupling prescription ([!]) to 



where (T) are the generators of SO(d) in the vector representation, with 
a = 1,2, ..,d, and (f) a (f) a = 1. In what follows, we shall restrict ourselves to 
the case of d — 3. In that case, the minimal coupling given by @ corresponds 
to gauging the vector 577(2) subgroup of the chiral SU(2) x SU(2) ~ 0(4). 

Note that the above gauging @ for d = 3 differs from that employed 
in ||. This is easily seen by identifying the SU(2) group valued field U 
as U = (j) a r a ,U- 1 = <p a f a , with r a = {ia a ,l), f a = (-ia a ,l), and A { = 
—iA^a a . In components then, it turns out that (D,jU) a = Di<j) a — Af<f> A and 
(DiU) 4 = D i( f) 4 + Af(f) a , in the notation of ©. 



2 The topological Invariant and the Model 

The cornerstone in our construction of the gauged Skyrme model on R3 is 
the selection of a suitable topological charge density which could be exploited 
to give a lower bound on the energy density of a suitably modified version 
of the 3 dimensional Skyrme model. Since we are dealing with a Skyrme- 
Sigma model featuring a constrained field and a global 0(4) invariance, we 
expect that the relevant topological charge density will be only a locally total 
derivative quantity as was realised in ||. The topological charge density of 
the 0(4) Skyrme H model is given by 





,d+l 



(2) 
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£o = e m e abcd d i( }> a di<l> b d k <t> c <l> d (3) 

where the index a = a, 4, and a = 1, 2, 3, etc., and the volume integral of (Q) 
yields the winding number provided that the fields exhibit the appropriate 
asymptotic behaviour. 

The gauged version of the density (^) , namely 

Qi = £ijk£ abcd Di(J) a Dj(f) b D k (f) c (f) d . (4) 

is given in terms of the d = 3 version of the covariant derivative @ 

D4 a = d i( f> a + e a ^A^\ A0 4 = <V- (5) 

The two densities (|]) and @ are related as follows 

Qi = £?o + Ze l]k ^F«D k <p a + «9A (6) 

where F°- is the 50(3) curvature of the connection defined by the first mem- 
ber of (g), and the density is given by 

a = 3e ijfe Q (2^<9 fc 4 - 4 Ff fc ). (7) 

The latter is a gauge variant quantity, like the density £>, (§). As we shall see 
below, the topological charge density will turn out to be 

g = go + ditti = q 1 - 3e ijk ^F«D k <j) a , (8) 

which is manifestly gauge invariant as it should be. 

The volume integral of the topological charge density g of (|j) in fact 
reduces to the volume integral of the winding number density g of (|3|), 
provided that the surface integral of Qi vanishes. This will be verified below 
for the spherically symmetric field configurations. Here we procede to explain 
our criteria by means of stating the following Belavin inequality: 

(AxA0° - ^e m e M D^D k ^f > (9) 
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where the square bracketsf..] on the indices imply antisymmetrisation, the 
constant K2 has the dimension of length and the constant A is dimensionless. 
The inequality is the first member of the two pairs, for the second member 
of which there are two options. These are, respectively, 



{k F<* + \g^)e ljk D k <\> a f > (10) 
(K Q h^)F« + L^D^f > 0, (11) 



in both of which the constant k q has the dimension of length and the yet 
undetermined functions g and h depend only on 4 , which according to the 
second member of (^) is the gauge invariant component of 4> a . 

As usual, the energy density will be the sum of the square terms of (|]) 
and (|T0|) , and respectively (|ll|), when these are expanded. Choosing the 



functions g(0 4 ) and h(4> 4 ) appropriately in each case, results in the cross 
term coinciding with the topological charge density q defined by (|8|). The 
topological charge, which is the volume integral of the latter, will be equal to 
the volume integral of g Q since the surface integral of fli vanishes as will be 
verified explicitly in the case of the spherically symmetric field configuration 
below. It then follows that the lower bound on the energy is the degree of the 
map, namely the winding number of the Hedgehog given by the volume inte- 
gral of £0, provided that suitable asymptotic conditions are satisfied. These 
are stated as usual to be 



lim 4 = -1, Jim 4 = 1. (12) 

x|^0 |:r|— >oo 

The appropriate choices for the functions g and h turn out to be the same, 
and namely, 

g^) = h(^) = 3A X (^)0 4 
respectively, leading to the two alternative Hamiltonian densities given by 
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Wi = ^(^) 2 +^(^) 2 (0 4 ) 2 (A0 Q ) 2 +^(A0 a ) 2 +^^P^A# 6 ) 2 (13) 
and 

H 2 = 9A 2 k 2 (</> 4 ) 2 (^) 2 + i(A0 a ) 2 + ^(A0 a ) 2 + \^{D^D 3 ^f. (14) 



We notice that both Hamiltonian densities break the global 0(4) symmetry 
of the corresponding ungauged sigma model, namely the Skyrme model. This 
is manifested through the appearance of |0 a | 2 and \Di(f) a \ 2 instead of \(f) a \ 2 
and \Di4> a \ 2 respectively. This situation can easily be altered by adding 
suitable positive definite terms to each of ( |i~3[) and ([L4]) without invalidating 
the respective topological inequalities. In the case of (|13|) for example, the 

quantity to be added is times (0 /3 ) 2 (A0 a ) 2 + (<%0 4 ) 2 . Both (|IJ) and 
flTj) then take the form 



n = riliF"? + |(A0 a ) 2 + (15) 

in which the constants rjo and 772 have the dimensions of length, and the 
constant t\ is dimensionless. The Hamiltonian density (ff3) coincides with 
the system proposed by Fadde'ev in Ref. |]]. 

The volume integrals of (0), ( |I4"D and ( pT5| ) are bounded from below by the 
winding number according to 



J d 3 xH > J d 3 xg = J d 3 xo , (16) 

provided that the field configurations satisfy at least the asymptotic condi- 
tions (H). 

It is not possible to saturate the inequalites ( |T6| ) by saturating the in- 
equalities (|jD and (0) separately to minimise absolutely the energy corre- 
sponding to (|T3"D , and, (0) and ([□]) separately to minimise the energy of (|1~4]) 
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absolutely. By minimising absolutely we mean solving the system by some 
first order Bogomol'nyi equations, and in the case of ( |l5f) this is not possible 
even in principle. For the models given by ( |T3"| ) and while in princi- 
ple possible, the saturated versions of (Qj and fllCf), ( pTTj) respectively, are 
overdetermined as is the case also for the (ungauged) Skyrme model[(||. This 
contrasts with the corresponding situation for the 2 dimensional 0(3) model 
where both the ungauged ||T0| and the gauged|| models support self-dual so- 
lutions. In the following therefore, we are concerned only with solutions of 
the second order Euler- Lagrange equations of the models (|13|), ( |14|) and (|i~5|) , 
and not with the solutions of some first order Bogomol'nyi equations. 



3 Spherically Symmetric Fields 

The spherically symmetric fields are given by the Ansatz 

(fi a = Xi sin/(r), 4 = cos/(r) (17) 

= atr)-l 
r 

yielding the following field strengths 

^ ,„ a sin f _„ , „. , a sin f s „ 

A0 Q = —^-5? + (/' cos / -^-)z<x Q (19) 

a' ( a' a 2 — 1 \ A A . , 

= 7 £ Va - ( 7 - J SijfiXpXa- (20) 

In ([T7|) -(POD x is the unit position vector and /' = df /dr etc. 

As stated in the previous section, we are concerned exclusively with the 
second order Euler-Lagrange equations here, and not in first order Bogo- 
mol'nyi equations obtained by saturating @ and fllPl) , and respectively @ 
and QTI|). These are easily verified in this spherically symmetric configura- 
tion (|T7| ) and ( |I"8|) to be overdetermined. We define the static Hamiltonian 
density of the one dimensional subsystem obtained by substituting ( |I9"D and 



6 



0|) into (jg), (0) and ([15]) respectively, by 

| dptf [/, /„; a, a p ] = ^J d 3 xH (21) 

where the dimensionless radial variable p is defined in the two cases as p = — 
and p = respectively, and f p — +. The definition ( |21~D is made, up to an 
unimportant constant multiple in each case. 



Since the alternative models (|T3|) , ( |i~4"D and (|TJ) are qualitatively similar, 
we shall restrict ourselves in the following to the detailed asymptotic and nu- 
merical study of one of these models only. We find it more natural to prefer 



models (0) and fll4|) to (15) because the former satisfy the minimal topolog- 
ical inequality fllEp, albeit without saturating it, while model ( |T5] ) satisfies 
an inequality derived from (|16|) itself. Next, we eschew model ([14]) because 
of its unconventional Yang-Mills term. Thus we restrict our considerations 
below to the model given by (|T3"|). 



In terms of the dimensionless parameter A2 = (^)j the resulting one 
dimensional Hamiltonian density for the models given by (pT3|) is 



(n 2 — 1 "l 2 Q \2 \ 2 

Hi = 2[2aJ + Z-^L] + cos 2 /[p 2 / 2 cos 2 / + 2a 2 sin 2 /] 

' p z 2 H 

+ |[p 2 / 2 + 2a 2 sin 2 /] + 2A 2 a 2 sin 2 f[2fi + «^±) (22) 
z p 

We first check that the topological charge, namely the volume integral 
of q given by ([§D reduces to the usual winding number, namely the volume 
integral of qq for the spherically symmetric field configuration ( |TTD and (PH) 
when the appropriate asymptotic conditions for the function f(r) 

lim/(r) = 0, lim/(r) = vr (23) 

are satisfied. In that case the volume integral of go is guaranteed to be the 
unit topological charge of the Hedgehog. 
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Concerning the asymptotic behaviour required of the fuction a(r) in the 
region r < 1, this is determined by regularity at the origin, while in the 
region r>l there are several possibilities consistent with a power decay of 
the function a(r) at infinity, which is a necessary condition for finite energy 
solutions. These asymptotic values are a(oo) = ±1 and a(oo) = 0. Unlike 
in SU(2) Higgs theory [[[]]], the behaviour of the gauge field and hence of the 
function a(r) at infinity is not directly relevant to the topological stability 
of the soliton. In the latter casejn| the topological charge coincides with 



the magnetic flux of the monoplole field, while here the topological charge is 
the degree of the map, which is the unit winding number for the spherically 
symmetric fields under consideration. Using the magnetic flux density Bi = 
^eijk<p a F^ k , for the spherically symmetric field configuration, we calculate the 
surface integral for the magnetic flux 

$ = J B.dS = A7t[(a 2 - 1) sin/] \ r=oa . (24) 

Irrespective of the value of o(oo), the flux <E> vanishes by virtue of the second 
member of (f2~3[). The models at hand therefore do not describe magnetic 
monopoles. 

Before stating our asymptotic conditions for the function a(r), we note 
that for the topological charge density g given by (§) to reduce to the usual 
winding number density go > the surface integral of f2j must vanish. This 
is seen by calculating the relevant one dimensional integrand, namely the 
quantity r 2 Xifli, from its definition (0), 

r 2 xA = (a 2 - l)sin2/ (25) 

which on the infinite 2-sphere clearly vanishes irrespective of which of the 
asymptotic values a(oo) = ±1 or a(r) = holds, provided that the corre- 
sponding condition stated by the second member of (^) does hold. 

Anticipating the results of our numerical integration to be carried out 
below, we state the asymptotic values of a(r) as 

lim a(r) = 1, lim a(r) = 0. (26) 

r->0 r^oo v ' v ' 

The Euler-Lagrange equations for the model given by (|22|) , with respect 
to the arbitrary variations of the functions f(r) and a(r) are respectively, 
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A 2 (l + 9A 2 cos 4 f){p 2 f p ) P + 8A 2 sin 2 / (a 2 f p ) p 

-18A 2 A 2 sin / cos / [p 2 f 2 cos 2 / + a 2 (cos 2 / - sin 2 /)] 

rp 1 sin ^ f 

+ 8A 2 a 2 sin / cos / [f 2 p - 2 J ] - 2A 2 a 2 sin / cos / = (27) 

and 

«pp - j 2 {a 2 - 1) - a sin 2 / [^A 2 A 2 cos 2 / + \\{ + A 2 (/ 2 + ^EB!/)] = 0, (28) 

The asymptotic values in the r«l region are given by the first members 
of fl23|) and (p6|), and for the model (|T3|) we find the following behaviours 

f(p) =7 r + Ap + o(p 3 ) (29) 

a(p) = 1 + Bp 2 + o(p i ) J (30) 
which lead, as expected, to differentiable fields (|T9"|) and (0) at the origin. 

The asymptotic behaviours of the solution in the r> 1 region are also 
power decays like for the usual (ungauged) SkyrmionQ. (Exponential decay 
can be obtained by incorporating a suitable gauge invariant 0(4) breaking 
potential ||.) Since we shall integrate the field equations with the asymptotic 
conditions (p3|) and (p6|) , we give the corresponding asymptotic solutions in 
this region 

f(p) = ~, <P) = ~, (31) 
P P p 

in which (3 = \{ ^/(C 2 A?(A? + 9A|) - 3) - 1], and the constants C and D will 
not be computed. 

Having solved the relevant Euler-Lagrange equations ( p6f) and (p7"D in the 
asymptotic regions p << 1 and p >> 1, we precede to integrate them nu- 
merically, subject to the asymptotic conditions fll8l) and fll9l). The constants 
A and .B in (|2"§|) and ([JO]) are fixed by the numerical integration. 
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The numerical integrations have been performed for the values of the di- 
mensionless coupling constants Ai = A2 = 1.5, and Ai = A2 = 1.4. The values 
of the pair of constants {A, B} were fixed in each of these cases respectibely 
to be 

{A,B} = {-1.9606513707554, -5.6289634247230} 

{A, B} = {-1.8663921477326,-4.7278588179841} 

The profiles of the function f(p) are given in Figure 1, and the profiles of the 
function a(p) in Figure 2. The profiles of the energy densities pertaining to 
each of these solutions are plotted in Figure 3, and they correspond to the 
total energies E 1 = 29.924879981245 and E 2 = 27.463710758882 respectively. 



4 Discussion and Summary 



Before proceding to summarise our results and making some qualitative com- 
ments, we give a breif quantitative description of the Baryonic current that 
the topological charge employed above pertains to. The latter is the volume 
integral of the density g given by which we identify with the fourth, 
time-like, component of this current j M . Accordingly, the full Minkowskian 
vector current is 



f = e^ vp(7 e abcd [D u cf) a D p (f) b D (T (f) c (f) d - h?f a D v ^$\. (32) 

The curvature field strength in ( |32"D consists only of an SU(2) field, say 
F"P = e^FJv, with F°^ = 0. Accordingly, the second term in ( |3"2"fl can be 
re-expressed using 



E^e abcd Ff a D v <p a <p h = 2e^°\2F«D v r<f ~ d v {F^ a 4>% (33) 



The second term on the right hand side of (|33|) being a total divergence, its 
volume integral vanishes and hence can be neglected. It follows from ([33]) 
that the divergence of the current j M given by (^) is 
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(34) 



The right hand side of (|34D can be shown to be locally total divergence which 
means that its volume integral vanishes and hence can be ignored, leading 
to a conserved Baryonic current, d^j^ = 0. This is exactly what we expect, 
since the vector gauging (|5|) does not lead to the divergence of the current 
being equal to an anomaly. Our current (^) can be compared to that of 
Goldstone and Wilczekjn], where the 0(4) sigma model has been gauged 



in the usual way according to D h 



and contrasted with 



the corresponding current of D'Hoker and Farhi fig] where the Skyrme model 
featuring the SU (2) valied field U has been gauged with the (V — A) SU (2) 
field according to D^U 
of the Baryonic current equals the anomaly. 



dfJJ —iA^.aU . In the latter case[]T5|], the divergence 



We have unfortunately not succeeded to adapt the constructions employed 
in the 0(4) sigma model with field <fr a , to the analogous case where the SU{2) 
valued Skyrme field U is used, which is the physically more interesting case 
as it leads to a non-conserved Baryonic current featuring the anomaly. Tech- 
nically, this has come about because of our inability to reduce the topological 
charge density ^ in this case, analogous to g used in the above, to the form 



q = so + d&i (35) 

in which <j = e^kTrll^diUdjU^dkU is equal to f? defined by (|3|) with 
U = <ft a T a , U^ 1 = (f) a f a defined at the end of Section 1. Had it turned 
out possible to establish (|3~5]) , then the topological charge would have been 
related to ft, analogous to Q\ in (||), which would have been the time-like 
component of the topological Baryonic current 



f = e^ pa Tr[U~ 1 D u U V~ x D p V \J- x D a U + -U^F^D^U]. (36) 



used in JT5[, whose divergence d^j^ does not vanish but is equal to the chiral 
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SU± anomaly. It would be very interesting if some other version of (|35| ) could 
be found, which would lead to a lower bound on the static Hamiltonian. 



To summarise, we have constructed three 50(3) gauged versions of the 
0(4) sigma model on R3, characterised by the static Hamiltonians (|13"D , (14]) 
and flT5|). These are equivalent to the corresponding gauged versions of the 
Skyrme model|4[]. Models (|T3| ) and (|l4])have the feature of breaking the global 
0(4) symmetry of the sigma model, but the latter can be restored by adding 
suitable positive definite terms to the Hamiltonian densities resulting in the 
model (p~5| ) which was first proposed by Fadde'ev|l[]. The large r asymptotic 
field configuration of the soliton presented here is S*0(3) symmetric and the 
magnetic flux of the corresponding field configuration vanishes. Accordingly 
this soliton model differs from the corresponding SU(2) Higgs model|]Tl|, in 
which the asymptotic field configuration is S0(2) symmetric and exhibits 
magnetic monopole flux. 



These models admit finite energy topologically stable soliton solutions, 
whose energy is bounded from below by the winding number, which can 
be interpreted as the Baryon number. This topological bound is saturated 
by Bogomol'nyi equations, which however are overdetermined as in the case 
of the usual (ungauged) Skyrme model, and hence likewise our solitons are 
solutions to the full second order Euler-Lagrange equations. These were 
solved analytically only in the asymptotic regions r << 1 and r >> 1, and 
the full integrations were performed numerically. 



In the present work, we have restricted ourselves to the spherically sym- 
metric case. As such the soliton in question carries Baryon number 1. It 
would be interesting to find the Baryon number 2 axially symmetric solu- 
tions, analogous to the corresponding axially symmetric solutions of the (un- 
gauged) Skyrme model|12[. Furthermore, since we know|| that for sigma 
models in odd dimensional spaces there are spherically symmetric solitons 
of arbitrary Baryon number N, it would be interesting to study these in the 
present model. These higher degree field configurations are characterised by 
their asymptotic values, which for small r differ from (^) according to 

Hm/(r) = 0, hm/(r) = Nn. (37) 

It would be interesting to integrate the Euler-Lagrange equations with the 
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asymptotic conditions fl3~T|), say with N=2, and to see whether the energy of 
that soliton is greater than twice the energy of the N = 1 soliton, as is the 
case for the usual (ungauged) Skyrme model|TB|]. All these detailed questions 
are deferred to future investigations. 
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Figure Captions: 

Figure 1. Profiles of the function f(p) for A = Ai = A2 = 1.5 and A = 1.4. 
The higher curve pertains to A = 1.5. 

Figure 2. Profiles of the function a(p) for A = Ai = A2 = 1.5 and A = 1.4. 
The higher curve pertains to A = 1.4. 

Figure 3. Profiles of the energy densities corresponding to the solutions with 
A = Ai = A2 = 1.5 and A = 1.4. The higher curve pertains to A = 1.5. 
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